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Abstract 

Minimax single facility location problems in multidimensional space 
with Chebyshev distance are examined within the framework of idem- 
potent algebra. The aim of the study is twofold: first, to give a new 
algebraic solution to the location problems, and second, to extend the 
area of application of idempotent algebra. A new algebraic approach 
based on investigation of extremal properties of eigenvalues for irre- 
ducible matrices is developed to solve multidimensional problems that 
involve minimization of functionals defined on idempotent vector semi- 
modules. Furthermore, an unconstrained location problem is consid- 
ered and then represented in the idempotent algebra settings. A new 
algebraic solution is given that reduces the problem to evaluation of 
the eigenvalue and eigenvectors of an appropriate matrix. Finally, the 
solution is extended to solve a constrained location problem. 

Key- Words: single facility location problem, Chebyshev distance, 
idempotent semifield, eigenvalue, eigenvector 



1 Introduction 

Location problems pQ form one of the classical research domains in opti- 
mization that has its origin dating back to XVIIth century and classical 
works by P. Fermat, E. Torricelli, J. J. Sylvester, J. Steiner, and A. Weber. 
Over many years a large body of research on this topic contributed to the 
development in various areas including integer programming, combinatorial 
and graph optimization (see, e.g. [2J [3J HI El [U E] ) ■ 

Among other solution approaches to location problems are models and 
methods of idempotent algebra El El EH HU Q2] , which find expanding ap- 
plications in the analysis of actual problems in engineering, manufacturing, 
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information technology, and other fields. Expressed in terms of idempo- 
tent algebra, a range of problems that are nonlinear in the ordinary sense, 
become linear and so allow more simple analysis and solution techniques. 
Specifically, many classical problems in graph optimization and dynamic 
programming reduce to solving linear vector equations, finding eigenvalues 
and eigenvectors of matrices, and to similar computational procedures. 

A single facility one-dimensional location problem on a graph is examined 
in [13} [8] , where it is turned into a problem of minimizing a rational function 
in the idempotent algebra sense. However, the proposed solution deals with 
polynomial and rational functions of one variable, and becomes less or no 
applicable in the multidimensional case. 

In [144 US] j a multidimensional constrained location problem on a graph 
is reduced to minimization of a max-separable objective function that can be 
represented as a maximum of functions each depending only on one variable. 
An efficient computational procedure is proposed which, however, seems 
to have limited application only to location problems where the objective 
function appears to be max-separable. 

In this paper, we further develop the algebraic approach proposed in 

IT71 [T8] . We consider a multidimensional minimax single facility loca- 
tion problem with Chebyshev distance, and show how the problem can be 
solved based on new results in the spectral theory of matrices in idempotent 
algebra. The aim of the paper is twofold: first, to give a new algebraic solu- 
tion to the location problem, and second, to extend the area of application 
of idempotent algebra. 

The rest of the paper is as follows. We begin with an overview of prelim- 
inary definitions and results in idempotent algebra, including basic concepts 
of scalar and matrix algebra, and elements of the spectral theory of matrices. 
Furthermore, a new algebraic approach based on investigation of extremal 
properties of eigenvalues for irreducible matrices is developed to solve mul- 
tidimensional problems that involve minimization of functionals defined on 
idempotent vector semimodules. 

We examine an unconstrained minimax location problem and represent 
it in terms of idempotent algebra. A new solution is given that reduces the 
problem to evaluation of the eigenvalue and eigenvectors of an irreducible 
matrix. Finally, the solution is extended to solve a constrained location 
problem. 

2 Preliminary Results 

We start with a brief overview of definitions, notation and preliminary re- 
sults of idempotent algebra that underlie the solution approach developed 
in subsequent sections. Further details can be found in [TJ El QUI [Til E2]- 
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2.1 Idempotent Semifield 

Let X be a set with two operations, addition © and multiplication ©, 
and their respective neutral elements, zero and identity 1. We suppose 
that (X, (D, 1, ffi, (8>) is a commutative semiring where addition is idempotent 
and multiplication is invertible. Since the nonzero elements of the semiring 
form a group under multiplication, the semiring is usually referred to as 
idempotent semifield. 

The integer power is defined in the ordinary way. Let us put X + = 
X \ {0}. For any x £ X + and integer p > 0, we have x — 1, (D p = 0, and 

x p = x v ~ x ®x = x ©x p_1 , x ~ p = (x~ 1 ) p . 

We assume that the integer power can naturally be extended to the case 
of rational and real exponents. 

In what follows, we omit, as is customary, the multiplication sign (g>. 
The power notation is used in the sense of idempotent algebra. 

The idempotent addition allows one to define a relation of partial order 
< such that x < y if and only if x @y = y. From the definition it follows 
that 

x < x (By, y < x ®y, 

as well as that the addition and multiplication are both isotonic. Below the 
relation symbols and the operator min are thought of as referring to this 
partial order. 

It is easy to verify that the binomial identity now takes the form 

(x © y) a = x a (By a 

for all real a > 0. 

As an example, one can consider the idempotent semifield of real num- 
bers 

K m ax,+ = (RU {-oo}, -oo, 0, max, +). 

In the semifield R m ax,+ , there are the null and identity elements defined 
as = — oo and 1 = 0. For each x £ R, there exists its inverse x~ l equal 
to the opposite number —x in conventional arithmetic. For any x,y € R, 
the power x y corresponds to the arithmetic product xy. The partial order 
induced by the idempotent addition coincides with the natural linear order 
defined on R . 

2.2 Vectors and Matrices 

Vector and matrix operations are routinely introduced based on the scalar 
addition and multiplication defined on X . Consider the Cartesian power X n 
with its elements represented as column vectors. For any two vectors x = 
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(xi) and y = (yi), and a scalar c G X, vector addition and multiplication 
by scalars follow the rules 

{x © y}j = Xj © yj, {ca;}j = cxj. 

A vector with all zero elements is referred to as zero vector and denoted 
by 0. 

The set X n with these operations is a vector semimodule over the idem- 
potent semifield X. 

A vector y G X n is linearly dependent on vectors x±,. . . ,x m G X ra , if 
there are scalars c\ , . . . , c m G X such that 

y = c\X\ © • • • © c m x m . 

In particular, the vector y is collinear with x, if y = ex. The zero 
vector is dependent on any vector. 

For any column vector x = (xi) G X™ , we define a row vector x = 
(x~) with its elements x^ = x~ l . For all x,y G X", the componentwise 
inequality x < y implies x~ > y . 

For conforming matrices A = (a^), B = (bij), and C = (cij) , matrix 
addition and multiplication together with multiplication by a scalar c G X 
are performed according to the formulas 

{A © B}ij = ciij © bij, {BC}ij = bikCkj, 

k 

{Cj4}jjf = CClij . 

A matrix with all zero entries is a zero matrix which is denoted by . 

Consider the set of square matrices X raxn . The matrix that has all 
diagonal entries equal to 1 and off-diagonal entries equal to is an identity 
matrix denoted by /. 

With respect to matrix addition and multiplication, the set X nxn forms 
idempotent semiring with identity. 

For any matrix A ^ and an integer p > 0, the power notation is 
routinely defined as 

A° = I, A p = A p ~ 1 A = AAP- 1 . 
The trace of the matrix A = (cnj) is calculated as 

n 

tr A = a i i . 
i=i 

A matrix is called irreducible if and only if it cannot be put in a block tri- 
angular form by simultaneous permutations of rows and columns. Otherwise 
the matrix is reducible. 
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2.3 Eigenvalues and Eigenvectors 



A scalar A £ X is eigenvalue of a matrix A € X nxn if there exists a nonzero 
vector x E X n such that 

Ac = Aa;. 

Any vector x ^ that satisfies the above equality is an eigenvector of A, 
corresponding to A . 

If the matrix A is irreducible, then it has only one eigenvalue given by 

n 

X= ^ti l / m {A m ). (1) 

m=l 

The corresponding eigenvectors of A have no zero entries and are found 
as follows. First we evaluate the matrix 

a x =A- 1 Ae-e(A" 1 A) n . 

Let a? be column i in A x , and be its diagonal element. Now each 
column af is replaced with that defined as 



a 



<, if 0^=1, 
0, otherwise. 



Furthermore, the set of columns is reduced by removing those columns, 
if any, that are linearly dependent on others. Finally, the rest columns are 
put together to form a matrix A + . 

The set of all eigenvectors of A corresponding to A (together with zero 
vector) coincides with the linear span of the columns of A + , whereas each 
vector takes the form 

x = A + v, 

where v is a nonzero vector of appropriate size. 



3 Extremal Property of Eigenvalues 

Suppose A G X nxn is an irreducible matrix with an eigenvalue A. For each 
x € X™ consider a function 

tp(x) = x Ax. 

It has been shown in |19[ [20] that (p(x) has a minimum equal to A and 
attained at any eigenvector of A. 

Now we improve this result by extending the set of vectors that provide 
the minimum of <£>{x) . 

First we revise the above result as follows. 
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Lemma 1. Let A = (a^) G X nxn 6e an irreducible matrix with an eigen- 
value X. Suppose u = (ui) and v = {vA are eigenvectors of the respective 
matrices A and A T . Then it holds that 

min x Ax = A, (2) 

where the minimum is attained at u and {v~) T . 

Proof. It is easy to verify that any vector x with nonzero elements satisfies 
the inequality x~ Ax > X. Indeed, let us take the eigenvector u and note 
that xu > (a^it) -1 / . Furthermore, we have 

x Ax = x Axu u > x^Au(x^u)^ 1 = X. 

It remains to present particular vectors x that turn the inequality into 
an equality. With x = u we have 

x Ax = u Au = Xu~u = X. 

Similarly, when x = (v~) T , we get the equality 

x~Ax = x T A T (x') T = v~A T v = Xv~v = X, 

which completes the proof. □ 

Assuming that a matrix A = (a^) € X™ x ™ is irreducible and has an 
eigenvalue A, we denote the set of vectors x that give minimum of x~Ax = 
A by 

Xa = arg min x Ax. 

Now we show that the set Xa is closed under main operations on vectors 
in X ra . 

Lemma 2. Suppose that x,y G Xa and c € X. Then the following state- 
ments are valid: 

(a) cx G Xa; 

(b) x®y(£ X A ; 

(c) {x~ ®y~)~ G X A . 

Proof. The first statement is obvious. To verify the next one, we take vectors 
x = (xi) and y = (t/i) , and consider a vector z = (zi) defined as z = x(&y. 
With the condition x, y G Xa , we have 

n n n n 

X = x-Ax © y Ay = 00 ^a^Xj © y^a^yj 

i=l j=l i=l j=l 

n n n n 

— 0^ ® 1(1 ij( X j © Vj) = Z i 1(1 ij Z j = Z 

i=l j=l i=l j=l 
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Hence we arrive at the inequality A > z~ Az . Since the opposite inequal- 
ity is always valid, we conclude that A = z~ Az . Therefore, z = x(By € Xa ■ 

The last statement is verified in much the same way. We put z = 
{x~ © y~)~ and then note that 



A = x Ax®y Ay = x i 1 i 

i=i j=i 



00% la ijVj 
i=i j=i 



> ©K rl © vi 1 )^ 1 ® yj 1 )' 1 = 00 *7 x «ij*i = z ~ Az - 
i=i j=i i=i j=i 

The rest of the proof is as before. □ 

Note that with the first and second statements of Lemma [21 the set Xa 
appears to be a vector subsemimodule in the semimodule X . 

Lemma 3. Suppose vectors x = (x^) and y = (yj) satisfy the condition 
x, y £ Xa ■ Then for all real a such that < a < 1 , it holds that 



l-a 



x lVl 



^a.X—a 
•^nan 



eX A . 



Proof. Assuming that z = (x^y\ a ,...,x"yj l a ) T , where < a < 1, we 
have 



A = (x-Axny-Ay) 1 - = 00 x~ a a^ y k 



- {1 ' a) alT a yl 



n n 

20© 



i=\ j=l 
-(l-a) 



k=l 1=1 

n n 



aijXj y- 



i=\ j=l 



1 ~ a - <^D £D 
j ~ WW 

i=l 3=1 



Using the same arguments as in the previous lemma, we arrive at the 
desired result z £ Xa ■ □ 

Consider a particular matrix A that has the form 

/ ai2 ... ai n \ 
a 2 i ... 



A 



\ a n i 



(3) 



/ 



where all entries au, ■ ■ ■ , «in and 021, . . . , a n i are assumed to be nonzero. 
Denote the first column and row of A as follows: 



a 



, Ol2, • • • j a ln) 
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Note that now we can write 

x~Ax = xa © b~x. 

It is not difficult to see that matrix is irreducible. In the case of this 
matrix, the result of Lemma can be refined as follows. 

Lemma 4. Suppose vectors x = (xi) and y = (jji) satisfy the condition 
x,y € Xa for matrix ([3]). Then for all real a« such that < a» < 1, it 
holds that 

: \(£X A . 

■^n tin / 

Proof. For each i = 1, . . . ,n we take a number a« such that < a» < 1, 
and define vectors 





( 


x 1 1 y 1 \ 




( x°t % y\ a% \ 


z = 






, Zi = 






{ 


rnOtnnX — Otn 

x n »n / 







It follows from Lemma [3] that Z\, . . . , z n € Xa ■ Furthermore, we write 

n n 

A = z i Az i = ©(*r a © h ~ z -) 



i=l 



i=l 
n n 



t=i j=i 



•///' <ij h t l xfy) - 



n 



z = z Az. 



i=l 



The rest of the proof goes through as before. 



□ 



By combining the results of Lemmas [JJ and [U we can arrive at the fol- 
lowing statement. 

Lemma 5. Let A be a matrix defined as ([3]) with an eigenvalue A. Suppose 
u = (v,i) and v = (uj) are eigenvectors of the respective matrices A and 
A T . Then ([2]) is valid for any vector 



x 



ai an— 1 

u 1 v 1 



< ai,...,a n < 1. 



Proof. It follows from Lemma[TJ that u G and (v ) T G X^. It remains 
to apply Lemma H] so as to complete the proof. □ 
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4 Unconstrained Location Problem 



In this section we examine a minimax single facility location problem with 
Chebyshev distance when no constraints are imposed on the feasible location 
area. 

Consider any two vectors r = (n, . . . , r n ) T and s = (si, . . . , s n ) T in R n . 
The Chebyshev distance (Loo or maximum metric) is calculated as 

p{r,s) = max \n — Sj|. (4) 

l<i<n 

Given m > 2 vectors rj = (rij, . . . , r n i) T € R n and constants iOj € K, 
i = 1, ... ,77t, the location problem under examination is to determine the 
vectors x = (x\, . . . , x n ) T G R n that provide the minimum 

min max (p(ri : x) + Wi). (5) 

a;eR" l<i<m 

Note that such problems are known as unweighted Rawls problems with 
addends [3J. Following the terminology of [2], the problem can also be 
referred to as the multidimensional Chebyshev Messenger Boy Problem. 

It is not difficult to solve the problem on the plane by using geometric 
arguments (see, eg, [HE]). Below we give a new algebraic solution that is 
based on representation of the problem in terms of the idempotent semifield 
fR m ax,+ , and application of the result from the previous section. 

4.1 Algebraic Representation 

First we rewrite @) as follows 

p(r, s) = s r © r s. 
Denote the objective function in problem ([5]) by f(x) and write 

m 

i=i 

With the vectors 

p = win © • • • © w m r m , q = wir± © • • • © w m r~ , 

we have 

m 

tp(x) = (£) Wi(x~ri © r~x) = x~p © q x, 

i=i 

and then represent problem ([5]) as 

min ip(x). (6) 

x£R n 
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Furthermore, we introduce a vector 



/ yo \ 




\ Vn J 

and a matrix of order n + 1 




Since we now have 

ip(x) = x~p © q x = y Ay, 
problem © reduces to that of the form 

min y Ay. (7) 

Note that the vectors y S R n+1 that solve ([7]) do not always have an 
appropriate form to give a solution to Specifically, to be consistent to 
([6]) , the vector y must have the first element equal to 1 . 



4.2 Algebraic Solution 

Consider problem ([7|), and note that the matrix A has the form of © and 
it is irreducible. It follows from Lemma [5] that 

min y Ay = A, 

where A is the eigenvector of A, and the minimum is attained at a vector 
that is obtained from eigenvectors u = (nj) and v = (vi) of matrices A and 
A T . 

First we evaluate A. For all k = 1, 2, . . . we have 

a-=w-p)"(; %), ^(^-'(y p «_), 

and therefore, 

ix(A 2k - 1 ) = 0, tr(^ 2fc ) = (q-p) k . 
Finally, application of (pQ) gives 

n 

A = tr^ m (A m ) = (q-p) 1/2 . 

m=l 
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To get vectors that produce the minimum in ([7]), we need to derive the 
eigenvectors of the matrices A and A T . Note that A T is obtained from A 
by replacement of p with (q~) T and q~ with p T . Therefore, it will suffice 
to find the eigenvectors for A , and then turn them into those for A T by the 
above replacement. 

To obtain the eigenvectors of A, we consider the matrix 

A-U = q 



p v 

Since for any A; = 1,2,... it holds that 



2k _ f „-„\-\ ( Q P U 
pq~ 



{\- l Af k = (q-p)- 1 
we arrive at the matrix A x in the form 

a x = A- 1 ^e---e(A- 1 ^) n+1 = ( 



1 (q-p)- 1/2 q- 
(q~p)~ 1/2 p [q p) X pq 



It is not difficult to verify that in the matrix A x , any column that has 1 
on the diagonal is collinear with the first column. Indeed, suppose that the 
submatrix {q~p) l pq has a diagonal element equal to 1, say the element 
in its first column (that corresponds to the second column of A x ). In this 
case, we have q p = c[\ Pi ■ The matrix A x takes the form 



A" 



4 1/2 -1/2 _ \ / fl -1/2 -1/2 

i ffi Pi q \ _ / i Qi Pi 

1/2 -1/2 -1 _ — 1/2 -1/2 -1 

Qi Pi p QiPi pq I V ^i Pi p Pi p 



where the second column obviously proves to be collinear with the first one. 

Let us construct a matrix A + that includes such columns of A x that 
have the diagonal element equal to 1 and are independent on each other. 
Since all the columns with 1 on the diagonal are collinear with the first one, 
they can be omitted. 

With the matrix A + formed from the first column of A x , we finally 
represent any eigenvector of A as 

u = ( 1/2 ) s, s £ R. 

V (<? p) 1 v J 

By replacing p with (q~) T and q~ with p T , we get the eigenvectors of 

A T 

{q-p)-\q-) T ) *' teR - 
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Application of Lemma [5] gives a solution of ([7]) in the form 



y 



(q-p) 1/2 ~ ai (pis) ai (qit) 1 - 01 



s, t G R, < ao, ... ,a n <1. 



With the condition that the first element of y must be equal to 1, we 
have to ensure the equation 

s ao t a -l = t 

to be valid for all «o such that < ao < 1 • Since the only solution to the 
equation is s = t = 1 , we arrive at the solution of ([6|) given by 



x 



< ax,...,a n < 1. 



\l/2— a nri a„ n l— an 



4.3 Summary of Results 

We summarize the above results in the form of the following statements. 
Lemma 6. Suppose that p = (pi) and q = (qi) are vectors such that 
p = win © • • • © u! m r m , q = u>irf • • • w m r~ n . 
Then the minimum in problem (|6|) is given by 

A = (q-p) 1 / 2 , 

and it is attained at the vector 

/ A Pi qi 

X = '■ 

i \l-2a„ a„ l-a„ 
\ /x Pn 'in 

for all ai such that < ctj < 1 , i = 1, . . . ,n . 

With the usual notation, we can reformulate the statement of Lemma [6] 
as follows. 

Corollary 1. Suppose that for each i = 1, . . . , n 

Pi = max(rji +wi,...,r im + w m ), qi = min(r;i - wi, . . . ,r im - w m ). 
Then the minimum in ([5]) is given by 

A = max(pi -qi,...,p n - q n )/2, 
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and it is attained at the vector 



( aiOi-A) 



x 



/ (1 -«!)(<?! + A) 



+ 



\ a n (p n - A) 
for all a, such that < on < 1, i = 1, . 



\ (1 - a n )(g n + A) 
n . 



An illustration of the solution is demonstrated in Fig. HH3J We start 
with two examples in Fig. [1] that present solutions in the plane R 2 when 
Wi = for all i = 1, . . . , m. In both examples, the given points are shown 
with thick dots, whereas the solution set is shown with a thick line segment. 



X 2 
P2- 



q-2. 
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Qi Pi 




Pi 



Figure 1: Solutions in R 2 when all Wi = 0. 



In geometric terms, the solution is obtained as follows. Construct a 
minimal upright rectangle enclosing all given points. Then trace two lines 
that are oriented at a 45° angle to the horizontal axis and go through the 
lower left and upper right vertices of the rectangle. The solution is the inner 
segment that these lines are cut off from the line drawn across the rectangle 
through the center points of its long sides. 

Examples of solution in the space R 3 are given in Fig. [51 where the 
solution sets take the form of rectangles depicted by thick lines. 

Fig. illustrates the solution of a problem with arbitrary constants Wi . 
First we present solution to an auxiliary problem obtained from the initial 
problem by setting Wi = for all i (top picture) . To get solution in the case 
of nonzero constants Wi , we replace each given point with two new points. 
Furthermore, the minimal rectangle is built for the new points and then the 
solution is derived in the same way as above (bottom picture). 
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Figure 3: Solution to a problem with nonzero Wi . 

5 Constrained Location Problems 

Suppose that there is a set S G R n given to specify a feasible location area 
in problem ([5]) and consider the constrained problem 

min max (p(ri,x) + Wi). (8) 

x£S l<i<m 

Representation in terms of the semifield K m ax,+ leads to the problem 

m 

min (I) Wip(ri,x). (9) 
i=l 

To solve the last problem we put it in the form of (0) by including the 
area constraints into the objective function of a normalized unconstrained 
problem. 
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5.1 A Normalized Problem 

First, problem © is transformed into a normalized form to enable subse- 
quent accommodation of the constraints in a natural way. We introduce 
new notation with a subscript 

p = Wl n © • • • © w m r m , = wirf © • • • © w m r~, X = (q^p ) 1/2 , 

and then define a normalized objective function 

<p Q (x) = Xq 1 {x~p © q^x). 

Instead of problem ([6]) , we consider the problem 

min lpq(x). 

xeR n 

It follows from Lemma [6] that the normalized problem has its minimum 
equal to 1 = 0, whereas its solution set obviously coincides with that of © . 

5.2 Maximum Distance Constraints 

Suppose that there are constraints imposed on the maximum Chebyshev 
distance from the facility location point to each given points. The constraints 
determine the feasible location set in the form 

S = {x G R n \p(ri, x) < di,i = 1, . . . , m}. 

For each i = 1, . . . , m, the inequality 

x~ri © r~x = p(ri, x) < di 

can be rewritten in an equivalent form as 

d^x'ri © d^r^x < l. 

With the notation 

Pi = d\ T\ © • • • © d^Tm, q^ = d^ 1 r 1 © • • • © d n ^r m , 

all constraints are replaced with one inequality 

x~p\ © q^x < t. 

Furthermore, we introduce a function 

ipi{x) = x pi © q^x 

and note that ip\ (x) < t if and only if the maximum distance constraints 
are satisfied. 



15 



Finally, we put 

p = Aq x pq © Pi , q~ = Aq 1 Qq © <h , 

and define the objective function 

ip{x) = ipo(x) © ifii(x) = x~p © q x. 

Now we can replace problem (J9j) by an unconstrained problem that has 
the form of problem ([6]) where the objective function <p(x) is defined as 
above. It is clear that both problems give the same solution set provided 
that the solution of the unconstrained problem has nonempty intersection 
with the feasible set. At the same time, the new problem allows one to get 
approximate solutions in the case when the intersection is empty. 

Based on the results offered by Lemma [U we can give a solution to the 
problem under the maximum distance constraints in the following form. 

Lemma 7. Suppose that p = {pi) and q = {qi) are vectors such that 
P = A 1 po©Pi, Q~ = Aq 1 9q © Qi", 

where 

Pa = w x r 1 ®---®w m r m , q$ = wirf © • • • © w m r~, A = (<?q Po) 1/2 , 
and 

Pi = d\ ffi • • • © dm r mi Qi = d-y 1 r 1 © • • • © d rr ^r m . 
Then the minimum in problem © is given by 

A = (q- P ) 1/2 , 

and it is attained at the vector 

\ 

for all cti such that < ai < 1, i = l,...,n. 

Going back to the usual notation, we arrive at the following result. 
Corollary 2. Suppose that for each % = 1, . . . , n 

Pi = maxOoi - A ,pii), Qi = min(<? i + A , qu), 
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where 



Poi = max(rji +w\,. 
q 0i = min(rji - w 1} . 
pu = max(rji - di,. 
qu = min(ra + dt, . 



• i Tim ~\~ ^m) j 

> Tim ^m) ; 

> Ti m d m ) , 
: Tim ~\~ d m ) , 



and 



,POn 



?0n)/2. 



A = max(p i - qoi, 
Then the minimum in ([9]) is given by 

A = max(pi - gi, . . . ,p n - q n )/2, 

and it is attained at the vector 
( ai(pi-A) 



x 



/ (!-«!)(</! + A) 



+ 



\ a n (p n - A) 



V (l-a„)(g„ + A) 



for all a.i such that < ati < 1, i = 1, . . . ,n. 

Fig. [5] gives an example of solution to a problem with maximum distance 
constraints in R 2 . The entire thick line segment represents the solution of 
the corresponding unconstrained problem, whereas the part of the segment 
inside the inner rectangle indicates the solution of the constrained problem. 
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Figure 4: Solution to a constrained problem. 
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